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Abstract. 

We report early BVR C I C imaging of the INTEGRAL GRB 041218. The observation started 129 s after the burst and a sequence 
of exposures in V, R c , I c , and B bands were repeated. The optical afterglow is detected in R c and / c bands while it was not 
detected in B and V. We find that the early afterglow deviates from a single power law decay and identify two characteristic 
timescales (breaks). Such breaks are theoretically expected from synchrotron cooling in the relativistic fireball (Sari et al. 1998) 
and we crudely constrain the physical condition of the relativistic shock. 
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1. Introduction 

The discovery of afterglow emission of gamma-ray bursts (van 
Paradijs et al. 1997 , Costa et al. 2004) opened new diagnostics 
for quantitative characterization of relativistic shock in GRBs. 
Follow-up observation within minutes of the gamma-ray emis- 
sion is particularly important, since the afterglow has informa- 
tion close (r < 0.01 pc) to the progenitor. Early optical after- 
glow has been detected only in a small number of occasions, 
including G RB 99 0123 (Ake rlof et al. [T9991 GRB 021004 
(Torii et al. 120021 Fox et al. 120031 . GRB 021211 (Wozniak 
et al. 120021 . and GRB 030418 (Rykoff et al. |2004i. In spite 
of the limited pieces of data, early afterglows have shown rich 
variety of behaviors which provided us with crucial informa- 
tion on the physical condition in the vicinity of the progenitor. 
These data brought us with the interesting pictures such as for- 
ward and reverse shock components (GRB 990123 and GRB 
021211), continued energy injection from the central engine 
(GRB 021004), and large extinction near the massive progeni- 
tor (GRB 030418). 

Since all these early detections were made by unfiltered or 
single filter CCDs, no information on the afterglow color could 
be derived in the previous studies. In this Letter, we report early 
BVR C I C imaging and the discovery of the optical afterglow of 
GRB 041218 with the Automated Response Telescope (ART). 
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2. Observations and Results 

A long (60 s) duration gamma-ray burst, GRB 041218 
was detected with the IBIS/ISGRI instruments onboard the 
INTEGRAL satellite at 2004 December 18, 15:45:25 UT 
(Mereghetti et al. 2004) and the burst position was distributed 
through IBAS (Mereghetti et al. 2003). After the earliest suc- 
cessful automated observation of GRB 030329 (Torii et al. 
120031 Vanderspek et al. 120041 . we upgraded the ART by using 
a larger telescope and installing standard photometric filters. 
GRB 041218 observations were performed with a 14 inch f/7 
Schmidt Cassegrain telescope, equipped with Cousin's BVRJ Q 
filters (Bessel 1990) in conventional filter wheel and cooled 
CCD camera (KAF-1001E). The first 60 s exposure in V band 
was started at 15:47:34 UT, or 129 s after the burst. Then a se- 
quence of 60 s exposures in R c , 7 C , and B bands was repeated. 
When the observation started, a large fraction of the sky was 
covered by clouds while the transparency gradually improved 
until 20 minutes post-burst. After that, the field of view was 
intermittently covered by clouds. 

We stacked four (two R c and two I c ) frames obtained on 
December 18 between 15:57:59 and 16:09:24 UT and visu- 
ally compared the resulting image with the digitized Palomar 
Observatory Sky Survey frame. In this stacked image, we noted 
a new object with S/N = 3.8 near the center of the 2'.5 radius 
error circle (Mereghetti et al. 2003 ) and reported it to the GCN 
Circular (Torii 2004 1. An independent identification of this new 
object and its fading was reported by Gorosabel et al. (2004). 
Figure □ shows the GRB 041218 field as imaged by ART and 
the digitized Palomar Observatory Sky Survey 2 Red for com- 
parison. 
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Table 1. Photometric results. 



Start Time 


Mean Epoch 


Magnitude* 


Filter 


Exposure 


(UT) 


(days) 






(s) 


15:47:34 


0.00184 


>12.19 


V 


60 


15:48:43 


0.00263 


>12.85 


R, 


60 


15:49:52 


0.00343 


>11.22 


1, 


60 


15:51:03 


0.00425 


>13.68 


B 


60 


15:52:12 


0.00505 


> 14.25 


V 


60 


15:53:21 


0.00585 


> 14.73 


Rc 


60 


15:54:29 


0.00664 


>15.33 


Ic 


60 


15:55:41 


0.00747 


>14.13 


B 


60 


15:56:50 


0.00827 


>15.29 


V 


60 


15:57:59 


0.00907 


16.25±0.54 


Rc 


60 


15:59:08 


0.00987 


15.43±0.27 


h 


60 


16:00:19 


0.01069 


>15.79 


B 


60 


16:01:28 


0.01149 


>16.18 


V 


60 


16:02:37 


0.01229 


>15.94 


R, 


60 


16:03:46 


0.01309 


> 14.72 


h 


60 


16:04:57 


0.01391 


>15.66 


B 


60 


16:06:06 


0.01471 


> 16.66 


V 


60 


16:07:15 


0.01550 


>16.68 


R, 


60 


16:08:24 


0.01630 


16.00±0.45 


h 


60 


16:09:35 


0.01712 


>15.21 


B 


60 


16:10:46 


0.01795 


> 16.66 


V 


60 


16:11:56 


0.01876 


>12.99 


Rc 


60 


16:42:10 


0.05425 


>17.62 


B 


60x9 


16:43:19 


0.05505 


>17.58 


V 


60 x 10 


16:44:28 


0.05585 


>17.50 


Rc 


60 x 10 


16:45:37 


0.05664 


>17.68 




60x10 



* Upper limits are at 2<x level. 



After dark subtraction and fiat fielding, we applied robust 
source detection and PSF photometry procedures to the data. 
When the source counts exceeded the 2<x level we considered 
it as a detection. In the other cases such 2<x upper limit was de- 
rived. The afterglow was detected in one R c frame at 2.2cr and 
two / c frames at 4.0 and 2.5<x level. The derived magnitudes, 
calibrated by using the field photometry of Henden (2004i, are 
summarized in Tabled Figure[2]shows the BVR C I C light curves 
combining our data with those reported in the GCN Circulars. 

3. Discussion 

The afterglow light curve of GRB 041218 shows, roughly, a 
power law decay. Here, we try to model the light curve as syn- 
chrotron radiation from the relativistic fireball. 

Fitting the R c band light curve with the function Mag = 
-2.5 ■ logf fl + b (tin days) gives a Q = -1.49 + 0.10 and b = 
23.47 ± 0.23 with ^ 2 /dof = 32.7/16. The errors for a and b 
are 90% confidence values. A large x 1 value is mainly due to 
fluctuations at 0.1 5 t J 0.3 days while this function marginally 
overestimate the flux in the first two R c detections, suggesting 
that the flux decay rate steepened at around 0.04 - 0.08 days. If 
we fit the R c data at t > 0.05 days, a steeper index a% = -1 .54 + 
0.12and£ 2 = 23.58+0.26are obtained with x 2 /dof = 30.2/14. 

For the I c band, only three detections are available and 
they are well described with a single power law with «3 = 
-1.49 + 0.16 and b 3 = 22.89 + 0.58. This function is drawn 
as dotted black line in Figure |5] for t > 0.00987 days. The 
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Fig. 1. Left: 6' x 6' field around the optical afterglow (at the 
cross point of vertical and horizontal bars) as imaged by ART in 
7 C band (2x60 s). The circle shows the 2'.5 radius INTEGRAL 
error region (Mereghetti et al. 2004). North is up and East is to 
the left. The narrow elongated object in the middle of the image 
is a cosmic ray event. Right: The GRB 041218 field in the dig- 
itized Palomar Observatory Sky Survey (DSS2 Red) is shown 
for comparison. A small circle shows the afterglow position. 
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Fig. 2. BVR C I C afterglow light curve of GRB 041218. The hor- 
izontal axis shows time from GRB onset. Blue, green, red, and 
black marks show, B, V, R c , and I c measurements, respectively. 
In addition to our measurements (filled circles; arrows mark 2cr 
upper limits), those reported in GCN circulars (Yonetoku et al. 
2004 Klotz et al. 2004 Ferrero et al. 120041 Monfardini et al. 
120041 D'Avanzo et al. l2004l Halpern et al. l2004l Cenko et al. 
2004 Mirsa & Pandey 2004 1 are plotted with open circles (de- 
tections) or arrows (upper limits). See discussion in text for the 
black and red lines. 

decay index in 7 C band is fully consistent with that for the R c 
band. However, the upper limit at 0.00664 days is significantly 
(0.59 mag) below the extrapolation from this power law. This 
means that the flux decay with 03 = -1.49 started just before 
the earliest detection at 0.00987 days. By using the upper limit 
at 0.00664 days, the decay index must have been a > -0.18 at 
t < 0.00987 days and suddenly changed to 03 = -1.49 within 
~ 5 minutes. 

We note that the observed changes in the decay indices (at 
0.00664 < t < 0.00987 days in I c and at 0.04 < t < 0.08 
days in R c ) are theoretically expected for synchrotron radia- 
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tion from relativistic fireball (e.g., Sari et al. 1998). According 
to this model, the early optical afterglow is expected to show 
breaks. The first one, t m , corresponds to the transit in the op- 
tical band of the typical electrons synchrotron frequency for 
minimum Lorentz factor, y m . The second one, t c , is the break 
corresponding to the transit of the cooling frequency in the ob- 
serving band. Therefore, we interpret t m =s 0.00664 - 0.00987 
days and t c - 0.04 - 0.08 days. 

We may constrain the power law index of electron number 
distribution, p, from the R c band light curve at t > 0.05 days. 
From the relation a = (2 - 3p)/4 and a 2 = -1.54 + 0.12, we 
derive p = 2.72 + 0.16. Before the cooling break at t c , the decay 
index is expected to be a — 3(1 - p)/4-. Therefore, we draw a 
power law decay with ci\ = -1.29 that starts with our first R c 
detection at t — 0.00907 days. These two power law functions 
are shown in Figure[2]by red lines; they intersect at t c 7.0 X 
10~ 2 days. 

We then constrain t m so that the limit in I a at t = 0.00664 
days is consistent with the model function. In the early part of 
the afterglow light curve, the flux is expected to increase with 
a power law index a = +1/2 (Sari et al. H998i . This constraint 
leads to a rising light curve at t 5 8.7 x 10~ 3 days, as drawn 
in Figure^ The existence of t m was originally invoked to con- 
serve the total energy above y m and it is unclear if such a clear- 
cut lower boundary exists in the accelerated electrons distribu- 
tion. The fact that the early afterglow seems to have brightened 
in a short timescale (-minutes) suggests that there is a distinct 
lower limit in the electron energy distribution at y m . 

For the other parameters, we can not solve the degeneracy 
due to the limited pieces of information. Instead, we try to build 
some trial parameter sets in the framework of the Sari et al. 
(1998 1 model with adiabatic evolution. An arbitrary assump- 
tion of isotropic burst kinetic energy E = 1 .Ox 10 52 ergs and the 
fraction of energy density in the magnetic field eg — 3.0 x 10~ 2 
gives the astrophysically reasonable set of other parameters, the 
fractional energy density of accelerated electrons, e e = 0.053, 
and the density of the surrounding medium, n = 4.2 ctrT 3 . 

The above discussion applies for constant density around 
the progenitor. The environment of long duration gamma-ray 
bursts are likely surrounded by stellar wind from a massive pro- 
genitor. In the framework of the Chevalier & Li (2000) model, 
the early flux rise in GRB 041218 afterglow may be explained 
as due to the clearance of the synchrotron self absorption. In 
this case (Case B in Chevalier & Li (2000)), the early after- 
glow brightens as oc f 7 ^ 4 . But again, the model parameters can 
not be fully constrained and we can not distinguish whether the 
interstellar or the circumstellar model gives the better descrip- 
tion. 

By interpolating R c and 7 C measurements, we may derive 
the afterglow color as (R c - 7 C ), =* 1.0 ± 0.6 at t = 0.00907 
days and (R c - 7 C ) 2 - 0.6 ± 0.3 at t = 0.1791 days. These val- 
ues are within the range (R - I = 0.46 ± 0.18, Simon et al. 
2001 1 measured in other afterglows. Although the Galactic ex- 
tinction in the burst direction is estimated as E(B - V)=0.61 or 
A(i? c )-A(7c)=0.45 from dust emission (Schlegel et al. 1998 1, 
the optical star count suggests a lower value (A(V)<1.0 or 
A(i? c )-A(/ C ) < 0.22) (Dobashi et al. l2005l figure 12). The fact 
that the GRB 041218 afterglow has an R c - I c color similar to 



that of other unextincted afterglows suggests that the lower ex- 
tinction is likely correct. The V - R c > -0.27 at t = 0.01471 
days is not that strong constraint and consistent with those ob- 
served in other afterglows (Simon et al. 2001 ). 

4. Conclusions 

We presented early BVR C I C imaging results and the discovery 
of the optical afterglow of GRB 041218. The afterglow light 
curve was found to deviate from a single power law decay and 
two characteristic timescales, t m =2 8.7 x 10~ 3 days and t c 
7.0 x 10~ 2 days, were identified. We showed that the standard 
fireball model gives simple enough account for the observed 
data. 

In the early phase of GRB afterglow, a variety of unex- 
plored physical processes, such as interaction of e ± pairs with 
the swept-up materials are expected (Beloborodov 2005). We 
expect that further BVRJ C follow-up observations with the 
ART for the INTEGRAL, HETE-2, and Swift GRBs will give 
insights into such phenomena. 
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